Transversity form factors of the pion in chiral quark models 
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The transversity form factors of the pion, involving matrix elements of bilocal tensor currents, 
are evaluated in chiral quark models, both in the local Nambu-Jona-Lasinio with the Pauli-Villars 
rcgularization, as well as in nonlocal models involving momentum-dependent quark mass. After 
suitable QCD evolution the agreement with recent lattice calculations is very good, in accordance 
to the fact that the spontaneously broken chiral symmetry governs the dynamics of the pion. Meson 
dominance of form factors with expected meson masses also works properly, conforming to the 
parton-hadron duality in the considered process. 
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The transversity form factors (TFFs) of the pion pro- 
vide valuable insight into chirally-odd generalized parton 
distribution functions (GPDs) as well as into the nontriv- 
ial spin structure of the pion. These interesting quantities 
have been determined for the first time on the lattice [l| . 
Formally, the TFFs, denoted as B^ ni (t), are defined as 
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where P' and P are the momenta of the pion, P = 
\{P' + P), A = P' - P, and t = A 2 . The symbol AS 
denotes symmetrization in v, . . . , fi n -i, followed by anti- 
symmetrization in /i, v, with the additional prescription 
that the traces in all index pairs are subtracted. The 
dividing factor of m„- is introduced by convention in or- 
der to have dimensionless form factors [lj. The tensor 
operators are given by 
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where D = ^(D — D), with D denoting the QCD covari- 
ant derivative. As in [l|, we use the positively charged 
pion and the up-quark density for definiteness. 

The available full-QCD lattice results [l| are for B*q U 
and -Bjj" and for —t reaching 2.5 GeV 2 , with moderately 
low values of the pion mass, m n ~ 600 MeV. The cal- 
culation uses the same set of QCDSF/UKQCD N f = 2 
ensembles with improved Wilson fermions and the Wil- 
son gauge-action that were used in the determination of 
the pion charge form factor [2j. 
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Form factors are related via sum rules to the GPDs (for 
extensive reviews see, e.g., |3l4ll| and references therein). 
Experimentally, the GPDs of the pion constitute rather 
elusive quantities which appear in rare exclusive pro- 
cesses, such as the deeply virtual Compton scattering 
or the hard electro-production of mesons. The high-Q 2 
dependence of the transversity form factors has been ad- 
dressed recently [121 ] . however the comparison with the 
lattice was avoided. In the present paper we fill this gap 
and confront the lattice transversity form factors with the 
results of chiral quark models, where particular attention 
is paid to spontaneous chiral symmetry breaking and the 
Goldstone nature of the pion as a composite relativis- 
tic qq bound state. We apply the Nambu-Jona-Lasinio 
(NJL) model with the Pauli-Villars (PV) regularization, 
as well as nonlocal chiral quark models inspired by the 
nontrivial structure of the QCD vacuum [13|, llJ]- These 
models provide the results at the quark- model scale 15J . 
After the necessary (multiplicative) QCD evolution 15J . 
our model results are in a quite remarkable agreement 
with the lattice data. Lower values of the constituent 
quark mass, ~ 250 MeV, are preferred. We use the tech- 
niques described in detail in [jja, [l6| . 

Previously, chiral quark models have proved to cor- 
rectly describe numerous features related to the pion 
GPDs. The parton distribution functions (PDF) have 
been evaluated in the NJL model in Refs. |17Hl9l |. The 
extension to diagonal GPDs in the impact parameter 
space was carried out in [201 ] . Other analyses of the pionic 
GPDs and PDFs were performed in nonlocal chiral quark 
models 0^, in the NJL model [11 Epl HHU and 
light- front constituent quark models 3l|, |32| . The par- 
ton distribution amplitudes (PDAs), related to the GPD 
via a low-energy theorem [33], were evaluated in [33 - l41| 
(see [42[ for a brief review of analyses of PDA) . The grav- 
itational form factors were computed in [43| ■ Finally, the 
pion-ph oton transition distribution amplitudes (TDAs) 
143-1471 were obtained in Refs. 48-521. 




FIG. 1. (Color online) The one-quark-loop triangle diagram 
contribution to the form factors B^ ni (t). 
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FIG. 2. (Color online) The NJL model parameters, m and A, 
plotted as functions of the constituent quark mass, M. 



In chiral quark models at the leading- N c level the cal- 
culation of the form factors and GPDs proceeds according 
to the one-loop diagrams (Fig. Q]), as explained in detail 
in [15|, [l6[ . The one-quark- loop action of the model is 
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where M is the constituent, and m the current quark 
mass. We apply the NJL with the PV rcgularization in 
the twice-subtracted version of Refs. [4l], |53, |54|. Vari- 
ants of chiral quark models differ in the way of performing 
the necessary regularization of the quark loop diagrams, 
which may to some extent influence the physical results \j 
Unlike many other studies, where one could work close 
to the chiral limit of m = 0, in the present case we need 
to tackle a situation with moderately large pion masses. 
This is because the lattice results for the transversity 
form factors are provided for m„ = 600 MeV. For that 
reason we do the following. As usual, the three model 
parameters A, M, and m are traded for the constituent 
quark mass, M, /„. (the pion decay constant), and m v . 
We assume that A depends on M only, and not on m. 
Constraining f n = 93 MeV (the physical value) and using 



1 We use the prescription where M 2 in the loop function is replaced 
with M 2 + A 2 , and then the regularized observable is evaluated 

O(0)-O(A 2 ) + A 2 (iO(A 2 )/dA 2 . 
= M 2 // 2 is not regularized. 



according to the formula O r cg ; 
The prc-multiplying factor g 2 



the given value of m^ leaves us with one free parameter 
only, M. The result of this procedure, with m for the 
two values of m^ of interest, is displayed in Fig. [2j 

An explicit evaluation of the one-quark-loop diagram 
of Fig. PR carried out along the standard lines explained, 
e.g., in [15j |. yields the simple result (holding at the quark- 
model scale): 
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with g-n = M/ fa and N c = 3 denoting the number of col- 
ors. The variables a and j3 are the Feynman parameters. 
Before comparing the results to the lattice data we 
need to carry out the QCD evolution, as the transversity 
form factors, not corresponding to conserved quantities, 
evolve with the scale. The lattice data correspond to 
the scale of about Q = 2 GeV, while the quark model 
calculation corresponds to a much lower scale, 



IM, = 320 MeV. 
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A detailed discussion of the evolution issue is presented 
in [151I55J 
tiplicatively as follows: 



It turns out that B^ Q and B^ Q evolve mul- 
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with the anomalous dimensions 'jTn = tt H n — 8 (H n = 
Ylk=i V 5 )) wn i cn gives 7ti = | and 7T2 = 8. We use 
A) = fN c - |JV> and a{n) = 4 7 r/[/3 log(/i 2 /A 2 , CD )], 
with Aqcd = 226 MeV and N c = N f = 3. In particular, 
this gives 
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Note a stronger reduction for Btiq compared to Btw- 
In the chiral limit and at t = 
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In Fig. [3] we show the results from the NJL model, 
evolved to /i — 2 GeV, confronted with the lattice data 
scanned from Fig. 1 of [l|. We have used m^ — 600 MeV 
and selected M — 250 MeV, which optimizes the com- 
parison. As we see, the agreement is remarkable. 

We have investigated the dependence of the values of 
the form factors at t — on ra^, as studied in [lj. The 
result is displayed in Fig. 0J We note a fair agreement in 
the intermediate values of m n , with a somewhat different 
character of the bent model curves and the flat data. 
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FIG. 3. (Color online) The transversity form factors ob- 
tained in the NJL model (lines) for M = 250 MeV and 
m^ — 600 MeV, evolved to the lattice scale of 2 GeV and 
compared to the lattice data from Fig. 1 of 1] (points). 
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FIG. 5. (Color online) The transversity form factors in the 
HTV model (solid line) and in the instanton-motivated model 
(dashed line) . The data as in Fig. [3] 
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FIG. 4. (Color online) The values of the transversity form 
factors at t — obtained in the NJL model (lines) for M = 
250 MeV and evolved to the lattice scale of 2 GeV, plotted 
as functions of m^ and compared to the lattice data of Fig. 4 
of [1| (points). 



Note, however, that the model, designed to work not too 
far from the chiral limit may need not be accurate at very 
large values of m w . Also, the lattice data are extrapolated 
to t = with a formula different from the NJL model, 
which may introduce some additional uncertainty. 

We have also explored the nonlocal chiral quark mod- 
els which incorporate the nontrivial structure of the QCD 
vacuum. In order to calculate the one-quark-loop dia- 
gram of Fig. [T]we use the nonperturbative quark propa- 
gator S(k) = 1/$ — m(k )] and the quark-pion vertex 



K(k,q) = ^ 5 T a F(k 2 + ,k 2 _), 



(10) 



where p± = k ± q/2. The quantity m (A: 2 ) is the dy- 
namical quark mass normalized by m(0) = Mq, and the 
nonlocal vertex F (A; 2 , fc 2 .) is normalized by F (fc 2 , k 2 ) = 
m (k 2 ). In the present study the nonlocal model calcula- 
tions are performed in the chiral limit, which means that 
m (k 2 -)• oo) = 0. 



Further, we will consider two variants of the quark-pion 
vertex (fTUl) . 



F I (k 2 + ,k 2 _) = Jm(k 
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F RTV (k 2 + ,k 2 _) = - [m(kl)+m(k 2 _)]. (12) 



The form (|10[) is motivated by the instanton picture of 
the QCD vacuum [13], while (fj"2"j) . the Holdom-Terning- 
Verbeek (HTV) vertex, comes from the nonlocal chiral 
quark model of [14]]. For t — both models yield the 
normalization 
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where m'(u) — dm{u) / du. In the local limit, where 
m(k 2 ) — > const, one reproduces Eqs. (|8I9[) . 

The results for B^ (t), n = 1, 2, are shown in Fig. [S] 
In the present study we have assumed that BT n o/ m n de- 
pends weakly on to„-, similarly to the local model (see 
Fig. S]). Hence, in order to compare to the lattice data 
for Bxno we simply multiply the results of calculations 
obtained in the chiral limit with m n — 600 MeV. We 
have carried out the same QCD evolution procedure in 
the nonlocal models as given by Eq. $3§. From Fig. [5] we 
note that the HTV model with the vertex function given 
by Eq. ([T2J (solid lines) and with M = 300 MeV works 
best, describing accurately the data, while the instanton 
model, Eq. (TTlT) (dashed lines), results in too steeply de- 
creasing form factors. Also, we have found that lower 
values of Mq spoil the agreement with the data. 
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FIG. 6. (Color online) Monopole fits to the transversity form 
factors. The bands correspond to the uncertainties of the 
parameters of Eq. (|16|) . The data as in Fig. [3] 



The ratio B^ (0)/BJ' 1 " (0) = A 2 /A x = 0.20(4) cor- 
responds, according to Eq. ©, to the evolution ratio 
a(fi)/a(iio) = 0.2(1), and hence to p = 350(80)MeVJn 
full agreement with the value ([5]) based on the PDF [17| 
and PDA [H of the pion (see 0, HI). 

The form factor BJ 10 couples to I G {J PC ) = 1 + (1~~ ) 
states, while B^ w to + (2 ++ ) and 1+(1 __ ) states. From 
Eq. (I16|) we note that indeed m\ is compatible with the 
mass of p(770), while rri2 with the mass of /2(1270), 
and within two standard deviations also with p(770) or 
p'(1435). These contributions cannot be disentangled 
with the current lattice accuracy. We note that the n = 2 
case allows also the coupling to the 1 + (1 H ) state, such 
as 6i(1235VjWhich, however, cannot decay into two pions 
(see, e.g., |59| for a discussion within Chiral Perturbation 
Theory). 
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FIG. 7. (Color online) Comparison of the predictions of the 
NJL model (solid line), the HTV model [3 (dashed line), 
and the monopole fit (dotted line). The data as in Fig. [3] 



In the large-iV c expansion all form factors are dom- 
inated by mesons with the proper quantum numbers 
(see, e.g., 56]). The well-known example is the exper- 
imentally measurable charge form factor, coupling to 
p(770), p'(1435), etc. (see, e.g., j57|), however meson 
dominance has also been checked in more elusive ob- 
jects such as the spin-2 gravitational form factor |43j 
(coupling to /2(1270)) and the trace-anomaly form fac- 
tor [58| (coupling to /o(600)). We thus undertake a sim- 
ple monopole x 2_n t to the TFF lattice data of 1] for 
600 MeV, reading 



•^TrcoW at "^ 



^TnOW — Ar, 



(15) 



and obtain 



A 1 = 0.97(6), mi 
A 2 = 0.20(3), m 2 



760(50) MeV, 
1120(250) MeV. (16) 



We conclude by presenting a comparison of the sev- 
eral considered chiral quark models in Fig. [7] We note 
the close proximity of all these model predictions. As 
we have shown, it is possible to describe the transversity 
form factors of the pion in chiral quark models. This is 
another manifestation of the fact that the spontaneously 
broken chiral symmetry is a key dynamical factor in the 
pion structure. Alternatively, one can describe the data 
with meson dominance, featuring parton-hadron duality 
for the TFFs. Indeed, appropriate meson masses gov- 
ern the fall-off of form factors, an expectation which be- 
comes exact in the large iV c -limit. The considered form 
factors, being the matrix elements of nonconserved cur- 
rents, undergo multiplicative QCD renormalization, thus 
their momentum dependence does not change as a func- 
tion of the scale, although the absolute normalization is 
governed by anomalous dimensions and the correspond- 
ing evolution ratio from the actual scale to the model 
reference scale. Actually, we find that the ratio of the 
lowest transversity form factors at t = is properly de- 
scribed when the QCD evolution is considered and the 
required model reference scale is fully compatible with 
other determinations. 
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